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Abstract 

Analytic formulae are derived for the complete photon energy spectrum due to QED 
corrections to fermion pair production in case of a limited angular acceptance for the final 
state fermions. After a numerical integration over the energy of non-observed photons, 
this corresponds to typical experimental conditions at LEP/SLC. 



1 Introduction 

One of the main tasks of the e + e~ storage rings LEP and SLC is a verification of the standard 
electroweak theory with unpreceded precision. An important reaction in this context is fermion 
pair production: 

e + + e- ^( 1 ,Z)-^f + + f- + n 1 , f ^ e. (1) 



Within the standard theory, the corresponding cross sections can be predicted with a precision 
well below 1 % either with (semi-) analytic formulae 0' @ or Monte-Carlo simulation 0. For 
the analysis of experimental data, the use of (semi-)analytic formulae is preferred due to 
their fast performance and well understood accuracy. A disadvantage is the small flexibility 
concerning the choice of experimental cuts. 

* Address to be expected after 3 October 1990: 0-1615 Zeuthen/Brandenburg(FRG) 
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In a series of papers, we derived analytic formulae taking into account QED corrections to 
the total cross section and the integrated forward-backward asymmetry as well as to the dif- 
ferential cross-section (see B and references quoted therein). The demand of experimentalists 
for fast and elegant algorithms motivated us to derive in addition analytic formulae for the 
photon energy spectrum in presence of an angular acceptance cut for the final state fermions. 
Such formulae allow a description of the quite realistic situation of fermion pair production 
including photon bremsstrahlung (with soft photon exponentiation) with cuts on the photon 
energy and on the fermion scattering angle using only one numeric integration. These new 
formulae are the subject of the present article. 

We denote with 0r(— ci, c 2 ) and Afb{— ci, c 2 ) the total cross section and forward-backward 
asymmetry, both integrated over the angular range cosO G (— c±, c 2 ): 

er T (-ci, c 2 ) = o- T (c 2 ) - or(-ci), (2) 

A F b(-cx,c 2 ) = ctt(-ci,c 2 )~ 1 1(Tfb(c2) + 0> B (-ci)]. (3) 

The building blocks a^(c),A = T,FB have to be determined by an explicit integration over 
the angular distribution: 

<ta(c) = E ^e[rf A (T^' (s, s;m,n) R a A (c;m,n)], (4) 

m,n=0,l a=e,i,f 

R A (c;m,n)=d A 1 [ A dv R A (c^m 7 n) a l ^ S ' ;m ' n l ^ 
Jo a A (s, s; m, n) 

4 

d T = -, d F B = 1, s = (1 - v)s. (6) 

Here, cr A °(s, s'; m, n) are normalisation factors containing the Born resonance function and all 
the couplings. Indices m,n are for photon (m = 0) and Z boson (m = 1) exchange. The 
remaining (numerical) integration in ( |5|) is to be performed over the photon energy with 
cut-off A G (0, 1). The problem to be solved is the determination of the angular dependent 
functions, 

R A (c,v;m,n) = / dcosO r A (cos6,v;m,n). (7) 

^0 

The general notations have been taken over from 0, where also the functions r A (cos8, v; m, n) 
may be read off. 

The following chapters are devoted to an explicit representation of rather compact formulae 
for the QED corrections ( |7j) originating from initial (a = e) and final (a = /) state radiation 
and their interference (a = i). 
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2 Initial State Radiation 

The major part of radiative corrections originates from the initial state. Therefore, an inclu- 
sion of higher order contributions is necessary. An adequate accuracy for LEP/SLC physics 
may be achieved if one combines the exact to order 0(a) cross sections with a soft photon ex- 
ponentiation procedure B 0. We split the QED corrections into two parts containing soft and 
hard bremsstrahlung, correspondingly. For initial state radiation, they both are independent 
of the kind of exchanged vector boson: 



R e A (c, v; m, n) = d A (c) [1 + S{f3 e )]l3 e v^ + H e A (c, v) 



d T (c) = d^(c + c 3 /3) 



T& + -Ql 
4 7T 



7T 

~3 



d FB (c) 

2 1 \ 



-a 8 
H%(c, v) = -Q 2 e d A l h A ( C) v)/v - f -d A (c), 

7T V 



2a 



71 



Ql(L e -l), 



In 



rnf 



A — T, FB, 



Qe = -1. 



(9) 
(10) 

(11) 
(12) 



The second term in ( [TTJ) is the remnant of soft photon exponentiation and makes the hard 
photon correction infra-red finite. The hard radiator functions are different for <jt and A FB ■ 
They depend on the reduced photon energy v and the range of observed scattering angles c: 



h e 



4 



T [c,v) 



-J 2 



7+ 7- 



z 3 ci 



+ 



2c 



7-7+ 



+2v z(c+c_) r 2 z — r 4 + 



zr 2 [(L e - 1) - lnz] 
4 z 2 22 , 



-z 2 (l - c + c_) + v 2 c + c_(z + r 2 c + c. 



3 c + c_ 



• 2v 4 (c+c ) :i ( -r,. - 2: J 



1 



: r 4 



(13) 



~2 ln 7T + 4^2 \ ~ 1) - M[2r 2 ^ 
rf 2 717+ ( H 



i> 2 



-r 2 c+c_ + Ar 2 c+c- — (z + r 2 c+c_)] 



r 



+w 2 c+c„[42; - 4c + c_r 2 + w 2 ] l-4zr 2 c + c_(^-^- + ^ ± ). (14) 

J 7- 7+ 

The following abbreviations are used: 

r n = l + z n , v = l-z, (15) 

c± = -(l±c), 7+ = c+ + 2c_, 7„=c„ + zc+. (16) 

For an unrestricted angular acceptance, c — > 1, the above formulae agree with those derived 
earlier for the uncut convolution kernels a F @' @ an d b . 
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3 Initial-Final State Interference Radiation 

Under usually realised experimental conditions, the 0(a) initial- final interference brems- 
strahlung is nearly negligible in the Z resonance region B H . This is not true if a tight cut 
on the photon energy is applied or in case of scattering angles near cos 9 « 1 . After an angular 
integration, the second type of singularity is removed completely. Away from the peak, the 
interference contributions become comparable to the other non-leading radiative corrections. 

It is sufficient to present the diagonal interference corrections. The non-diagonal terms 
obey the following simple relation: 

R A (c,v;m,n) = ^[R A (c,v;m,m) + R A (c,v;n,n)*}, A = T, FB. (17) 

The initial-final state interference contributions are composed of soft and hard bremsstrahlung 
parts S A (c, e), h A (c, v) and of 77 and 7Z box diagrams Ba(c; m, n): 

R A (c,v;m,n) = -Q e Q f {5{v)[S A (c,e) + B A (c;m,n)] + 6(v - e)h\(c,v)a\(s, s';m,n)}, (18) 

7T 

where A = T,FB if A = FB,T. The box contributions are the only ones with an explicit 
dependence on the kind of exchanged vector boson; they are written in (anti-) symmetrised 
form: 

B t ,fb(c] m, n) = b TjFB (c; m, n) ± b TjFB (-c; m, n). (19) 
The corrections Rip(c,v; m,n) to the total cross section cr T (c) contain the following terms: 

S l T (c, e) = -4 lne[(c 2 - 1) ln^± + 2c] + 2(c 2 - l)[Li 2 (c+) - Li 2 (c_)] 

-(c 2 - l)ln(c+c_)ln— -4cln(c+c_) -4 In — + 8c, (20) 

c_ c_ 

b T (c; 0, 0) = |(c 2 - 1) In 2 c+ + 2clnc+ - (c 2 - 3) lnc+ - 3c - i[vr(c 2 - 1) lnc+], (21) 

b T (c;n,n) = -21n(l - -){2cR{l - R) lnc+ 

R 

-[-2R 2 + R{c 2 + 1) + c 2 - 1] lnc+ - cR{R + 1)} 

+AR{R - l)cLi 2 (l - h - 2c\nR+ (c 2 - 1) In 2 c+ 

R 

+4clnc+ + 2[R(c 2 - 1) - c 2 + 3] lnc + - AcR(R - \)Li 2 {\ - 

R 

-2[2 J R 2 - J R(c 2 + l)]Li 2 (l-^)+2 J Rc-6c, n ^ 0, (22) 

R 

R = - s [M 2 z -iM z T z (s)l (23) 



h l T (c, v) = 4c + c_ 



r -^\n C -±+(^z-v 2 )\n^± 

v c_ V c + c_ / 7_ 



+ 4cz (-— +lnz) . (24) 
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The interference corrections R % FB (c,v\m,n) to the numerator of the forward-backward 
asymmetry, (Jfb(c) , contain the following terms: 

^S FB (c, e) = In e[-8 In 2 - 4 ln(c+c_) - 3(y + c) In °± - c 2 ] 

+~(|- + c)[Li 2 {c + ) - L* 2 (c_)] + 2[Li 2 {c + ) + Li 2 {c^)} - ~(1 + c ) ln(c + c_) In ^ 

1 c 2 
-(In 2 c+ + In 2 c_) + 4 In 2 2 + In 2 - - ln(c+c_)(c 2 - 1) + 2Li 2 (l), (25) 

6 FB (c; 0, 0) = — (c 2 - 1) In 2 c+ + (c 2 - 3) In c+ - 2c In c+ - ^(ln 2 2 + 6 In 2 + c 2 ) 

-itc[^ In 2 + (c 2 + |) lnc + + 2(^ + c) lnc + - ^c 2 ], (26) 



6fb(c; n, n) = (— + c) In c+ + - In c+ 



+ ln(l - i){ (4i? 2 - 2i? + H) l n 2 + [4i? 2 - 2i?(c 2 + 1) + 2c 2 + H] l nc+ 
-fx I 3 3 

+4[cfl(.R - 1) + (y + c)] lnc+ + c 2 (-# 2 + 3i? — |)| 
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+c - l - B -_ji„ B+lifl 3_ 4B2 
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4 
3 



2 "> l)-^ln 2 2 + (^ 2 + ^-6)ln2 



+ 



+2c\R-l)Li 2 (1 
4 



-i? 2 + i?(— c 2 + -) + 2(c 2 - 3) 
3 V 3 3 J K J 



\nr + \ ( ^R 2 + -R-4) rhir_ 



+2 



--i? 3 + 2i? 2 -i?(c 2 + l) + c 2 
3 



+2 [2c(i? 2 - R) + ( °- +c)] U 2 (l - ^) + " 1)' 



Li 2 [ 1 - 

2 



<"■+ 



R 



(27) 



,« / N 2c2z f l 2 2 2 ^ 1 2 

ft FB (c, v) = I -« - — + — + -C Wi 

7+7- V3 3f ri 3 

— 2c 2 zu ln z — 4 — c + c_r 2 ln(c + c_) 
f 

2 



- — — (ci lnc_ + cl In c_|_) H — [ 2z — 5r 2 ) — lnri 
3 f 



-4 ln 7_ 
+4 In 7+ 



n 

v 

2CC 2 . -7_(wH Ac 2 _ +7_ 

3 f 

2cc + + 7+ (?; + -^± -4c 2 + 7+ 
3 w 



(2f 



[0 



Again, for an unrestricted angular acceptance the above formulae yield those derived earlier 
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4 Final State Radiation 



The final state radiator functions for the angular distribution, r A (cos6, v; m, n), as used in ( |?]) 
for the definition of R^(c,v;m,n) are described in detail in chapter 5.1 of ®. Their angular 
dependence to order 0(a) is the same as that of the Born cross section. Consequently, an 
integration is trivial. For applications which deserve highest precision, a common exponen- 
tiation of initial and final state radiation corrections is recommended (see also 1 0). An 
inspection of the formulae in chapter 5.3 of @ shows that the only nontrivial angular depen- 
dence is due to the initial state hard radiator function. The result of their integration may 
be found in section 2 of the present article. Further, there exists a minor angular variation in 
the non-leading part of the final state factor which is proportional to [3 — 4/(1 + cos 2 9)]. We 
have shown numerically that this term can be neglected completely for all applications. 

As a result, a description of common soft photon exponentiation of initial and final state 
bremsstrahlung is obtained by a direct combination of the treatment in with the integrated 
expressions for the initial state hard photon corrections ( ( [14]). 



5 Discussion 

The formulae of this article may be used both in a model-independent way with input pa- 
rameters chosen to be e.g. M z , T z , T e , Tf, or within the standard theory using e.g. a, GL, M z 
[0> 0. The corresponding Fortran codes ZBIZON and ZFITTER are described in and 



have been applied recently to data obtained by LEP experiments (e.g. 1^1 |20| y j n a 
unique way, one can perform a semi-analytical calculation of either observables which are in- 
tegrated over a wide angular range or of pseudo-differential distributions with bins filled using 
differences of angular slices. As typical examples, we show in Figs. 1 and 2 some predictions 
for such distributions of o~t and A FB . 

To summarise, we think that the analytic calculation of QED corrections with partial 
angular integration could prove to be a powerful tool for the phenomenological analysis of 
fermion pair production. 



1 Besides the calculational chain based on the formulae of this article, these codes contain also a branch 
with different choice of cuts ILij . 
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List of Figures 

Fig. 1. 

The total cross section a T {c\, c 2 )/Ac, Ac = c 2 — Ci, as function of the scattering angle and of 
the photon energy cut-off A = 2E 1 / y /s] M H = 100 GeV, m t = 100 GeV. 

Fig. 2. 

Cumulative forward-backward asymmetry Afb{— c, c) as a function of c and of the photon 
energy cut-off; parameters as in fig. 1. 
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